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Abstract 

Dark matter decaying or annihilating into /i + jJL~ or t + t~ has been proposed as an 
explanation for the anomalies reported by PAMELA and Fermi. Recent analyses 
show that IceCube, supplemented by DeepCore, will be able to significantly constrain 
the parameter space of decays to u + fi~, and rule out decays to t + t~ and annihila- 
tions to fi^ u~ in less than five years of running. These analyses rely on measuring 
track-like events in IceCube+DeepCore from down-going v^. In this paper we show 
that by instead measuring cascade events, which are induced by all neutrino flavors, 
IceCube+DeepCore can rule out decays to u + u~ in only three years of running, and 
rule out decays to t + t~ and annihilation to u + u~ in only one year of running. These 
constraints are highly robust to the choice of dark matter halo profile and independent 
of dark matter-nucleon cross section. 



1 Introduction 



The existence of dark matter has been established by numerous observations. Although it 
constitutes most of the matter in the universe [1] , the nature of dark matter remains largely 
unknown. One widely held possibility is that it is a new fundamental particle produced in the 
early universe and present today as a thermal relic. Among the best motivated of these are 
the so-called weakly interacting massive particles, or "WIMPs" (for reviews, see Refs. [21 E]) 
which are predicted to be undergoing annihilations [U El El E] and possibly decays [8] in the 
current epoch. 

Recently, the instruments PAMELA [9J , ATIC [10J, PPB-BETS [H], and Fermi p| have 
observed features in the cosmic-ray e spectrum and a positron fraction that are inconsistent 
with known backgrounds. While these anomalies may be due to unidentified astrophysical 
sources [13], one exciting possibility is that they are due to the decay [H f [T5 ] \W \ [T7 1 [TBI \19\ 120] 
or annihilation [211 E21 E3J [2H [25j |26] of dark matter particles into standard model states. 
Even in the case that anomalies are not caused by new physics in the dark sector, the 
constraints are generally applicable to dark matter models. 

In order for dark matter to explain the anomalies, the products of decay or annihilation 
must be primarily leptonic. In either case, the Fermi observation gives the most precise 
preferred region in mass and lifetime/cross-section; for decays, only fi + fi~ and t + t~ final 
states fit the data, while for annihilations a small region of e + e~ is also permitted [29] . For 
the allowed parameter space, decays or annihilations into hadrons and weak and Higgs bosons 
are severely limited. Decays are constrained by the Fermi observation of the isotropic diffuse 
gamma-ray flux [31] , and annihilations are constrained by the production of antiprotons [32J . 

Explaining the anomalies with annihilations has an additional challenge. In order to 
match the observed rates, the cross-section required is 10 3 -10 4 times that expected for ther- 
mal production of the dark matter in the early universe. This necessitates nonthermal 
production mechanisms or low-velocity enhancements to the cross-section, such as the Som- 
merfeld enhancement [211 12S] or the Breit-Wigner enhancement [331 El] ■ Moreover, high cross 
section annihilations to leptonic states are tightly constrained by synchrotron radio emis- 
sions from the galactic center, although this constraint can be evaded if the true galactic dark 
matter halo profile is much less steep at the galactic center than benchmark profiles fI7\ |2"U] . 
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For either decays or annihilations, neutrino observations may provide the strongest con- 
straints — or the most promising corroborating signatures — for dark matter to be the source 
of the anomalies. With the exception of the disfavored e + e~ channel, the required leptonic 
final states decay into neutrinos. These travel undeflected from their sources, eliminating 
any uncertainties in modeling their propagation. Moreover, astrophysical sources that may 
explain the anomalies are not expected to produce a large flux of neutrinos. Pulsars, for 
example, would generate at most 0(1) events/year at a full km-scale detector, and only at 
energies greater than 10 TeV [35J. 

There are, however, challenges to constraining dark matter models with neutrinos. As 
they are observed only by collecting Cherenkov light from induced particle showers or from 
secondary muons, angular resolution is poor compared to gamma-ray observatories. Also, 
there are significant backgrounds from atmospheric muon and neutrino fluxes. However, 
if these backgrounds can be controlled, the poor angular resolution need not be a barrier; 
indeed by integrating the flux over a large area of the sky, the resulting constraint is much 
less sensitive to the choice of dark matter profile [39j HQ] . Moreover, by observing the galactic 
and extragalactic diffuse dark matter, rather than any that may have been captured by the 
Sun or the Earth, any constraints will be independent of the dark matter-nucleon cross 
sections, which are related to final states in a model-specific way |41j . 

Recent analyses show the power of neutrino constraints, using various strategies to re- 
duce the effect of atmospheric backgrounds. The Super-Kamiokande observatory resides 
in the northern hemisphere, facing away from the galactic center, minimizing atmospheric 
backgrounds. Measurements of upward-going muons place a limit on the flux of galactic u^, 
providing a robust constraint that eliminates annihilations to t + t~ as a source of the e 
anomalies [281 EH]. The IceCube observatory, on the other hand, resides at the South Pole 
where down-going atmospheric fluxes are coincident with the neutrinos from the galactic 
center. The overwhelming background of atmospheric muons can be suppressed by event 
selection to establish an isotropic diffuse flux limit [321 231 133] , but this limit only starts at 
a high energy threshold O(10 TeV), and yields a relatively weak constraint as we will show. 

Currently under construction is DeepCore [33], an in-fill of the IceCube detector which 
will use the outer instrumented volume as a veto on downward-going muons to a level of 
one part in 10 6 [46J . This will allow the galactic neutrino flux to be measured and compared 
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against the atmospheric neutrino flux, providing a constraint on dark matter decays and 
annihilations. Recent work [4"Tl |4*8] by some of the authors shows that IceCube+DeepCore 
will be able to significantly constrain the parameter space of decays to , and rule out 
decays to r + r _ and annihilation to as possible sources of the anomalies in less than 

five years of running. 

Recently, IceCube+DeepCore has demonstrated in simulations the ability to distinguish 
between track-like events, which are due to the charged-current interactions of u^, and cas- 
cade events, which are induced by v e ^ T through charged-current interactions and by all neu- 
trino flavors through neutral- current interactions [4*9ll50] . This is very useful for constraining 
dark matter neutrino fluxes because is the dominant flavor of atmospheric neutrinos above 
40 GeV (511 152] . The neutrino-nucleon cross sections are the same for all flavors, so signal 
would be considerably enhanced, while background would be reduced because only creates 
cascade events through the neutral-current interaction, which is lower in cross section than 
the charged-current interaction J37J EH E9] (see Fig. [I]). Moreover, cascade events are easy 
to distinguish from the tracks caused by any atmospheric muons that are not vetoed by the 
outer volume. 

In this paper we show that observation of cascade events at IceCube+DeepCore can 
enhance the neutrino constraints on dark matter, and rule out (or corroborate) dark matter 
decays to or t + t~ and annihilation to fi + as sources of the observed e cosmic-ray 

anomalies in a much shorter time compared to searches that rely solely on track-like events. 

2 Analysis 

For brevity, in the expressions below we write only the terms for neutrinos and not the terms 
for antineutrinos. The terms are identical, except replacing "z/" by "z/". 

2.1 Galactic flux signal and background 

The flux of neutrinos from the galactic dark matter halo is given by 
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Figure 1: Cross sections for (anti)neutrino-nucleon interactions, given by Ref. [37]. The 
blue lines are for charged- current interactions, and the red lines are for neutral-current 
interactions. Solid lines are for neutrinos and dashed lines are for antineutrinos. 

for annihilations, where r Q = 8.5 kpc is the distance from the Sun to the galactic center [53J, 
Pq = 0.3 GeV cm -3 is the dark matter density in the solar neighborhood, m is the dark 
matter mass, and r and (av) are the dark matter lifetime and thermally averaged cross 
section respectively. J n is the line-of-sight (los) integral through the halo profile: 

UAQ) = ^- f dSlf -(■£-)" (3) 
AiZ J An Ao.s. r Q \PeJ 

where Afl is the region of sky observed. In this analysis we use the Navarro-Frenk- White 
(NFW) halo profile (55] 

with r s = 20 kpc. 

The neutrino source spectra dN v JdE for flavors i are given by PYTHIA [56J simulation. 
Assuming sin 2 2#i2 : sin 2 26 l i3 : sin 2 2#23 ~ 1 : 1 : 0, the flavor distribution will be 1 : 1 : 1 as 
the neutrinos reach the Earth, having traveled a variety of long distances across the galaxy 
and being well mixed through vacuum oscillations. 

As discussed in the Introduction, the only background to cascade events is atmospheric 
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neutrinos after the veto suppresses the background of atmospheric muons to one part in 10 6 
and event selection is used to eliminate the rest. The and v e fluxes are given by Ref. [51], 
where the flux agrees well with AMANDA observation [54J . At low zenith angles, the flux 
of background is ~ 20 times greater than the flux of v e from 40 GeV to 100 TeV; it is 
~ 1000 times greater than the flux of v T (see Ref. [52]). However, below 40 GeV, especially 
at high zenith angles, the atmospheric fluxes of the three flavors only differ by 0(1) due 
to flavor mixing. Since modeling the background (and signal) below 40 GeV would require 
simulating flavor mixing as the neutrinos propagate through the limb of the Earth, for this 
analysis we simply set an energy cutoff of -^thresh = 40 GeV. 

To obtain the event rates due to the galactic and background fluxes, we first set energy 
bins for dark matter of mass m to be [max(i? t h re sh; m/10),m/2] for decays and [max(.E t h r esh, m/5),m] 
for annihilations. Note that the bin width is much greater than the energy resolution, 
log 10 (-E max /-E min ) ~ 0.4 for track-like events and log 10 (£ , max /£ , min ) ~ 0.18 for cascade events 
We then integrate the flux times the effective area over the energy for each bin. 

For track-like events due to v^, the event rate is 



-'max 



r tr . = I dQ I dE p icc N A V tI . [ [<r vN (E)]cc (5) 



where p- lce = 0.9 gem -3 is the density of ice, Na = 6.022 x 10 23 g -1 is Avogadro's number (to 
convert mass to number of nucleons), [o" i/ Ar(-E')]cc is the neutrino-nucleon cross section for the 
charged- current interaction, and V tI . ~ 0.04 km 3 is the effective volume of the detector for 
track-like events [46]. Note that we do not add the residual atmospheric muon background 
to the background of track-like events due to the uncertainty in its value. 

For cascade events we use the instrumented volume V casc ~ 0.02 km 3 [HUE?], the charged- 
current interaction for z/ 6)T , and the neutral-current interaction for all flavors to obtain 

Tease. = J dn j E dE P^N A V caac . [[<T vN (E)]cC + Kn(E)] NC J " ( 6 ) 

Unlike track-like events, cascade events are well contained, so the effective volume for their 
detection will vary little with energy. Also, we assume that in neutral-current interactions 
all the energy of the neutrino goes into the cascade. Taking partial energy transfer into 
account would yield a modest improvement in significance, since most of the signal is from 
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v e>T charged-current interactions but most of the background is from neutral-current 
interactions J3 

Finally, because the pointing capability for cascades is approximately 50° [46, 49J and 
the pointing capability for track-like events has yet to be established, we integrate over the 
2tt field of view around the galactic center. As mentioned before, this should provide a 
constraint which is robust to changes in dark matter halo profile. Specifically, the fractional 
change between the NFW profile used here and the much less steep isothermal profile for 
Ji(2tt) is 0(l(r 3 ), and for J 2 (2vr) is 0(0.1). 

2.2 Extragalactic isotropic diffuse flux 

For comparison to the DeepCore constraints from down-going fluxes, we calculate the con- 
straint from the isotropic diffuse flux limit for AMANDA-II from track-like events [43J and 
the projected limit for IceCube from cascade events |44j. 

The main contribution from dark matter decay to the isotropic diffuse flux is that from 
extragalactic dark matter residing on cosmological scales. The flux from the decay of cos- 
mological dark matter is given by the previously derived (301 El] formula 

( A (ex ° = c n DM p c i r y- 3/2 ( dN Vi \ 

\dEdQj " 4n Ho Q^mr J, V ^l + n K /Q M y^ \d{Ey) J 1 } 

where y = 1 + z, with z being the redshift, Hq = 71.9 kms -1 Mpc -1 is the Hubble constant, 
p c = 3Hq/(8ttGn) is the critical density, f2 M = 0.258, f2 DM = 0.214, and = 0.721 are, 
respectively, the total matter, dark matter, and dark energy densities divided by the critical 
density [lj. The isotropic diffuse flux from the annihilations of cosmological dark matter is 
too small to be relevant, since the density of dark matter on cosmological scales is very low, 
and the flux is suppressed by another power of p c /m. 

Because of the loss of signal due to event selection and the presence of background fluxes, 
both the AMANDA-II limit and the projected IceCube limit are only valid at energies greater 
than ~ 20 TeV. Below these thresholds we add the atmospheric background flux to these 
limits, and use these total fluxes to calculate the constraints on the dark matter lifetime. 



Inelasticity curves are given in Ref. but only for energies 10 TeV and higher. 
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Figure 2: Constraints for decay to [i + fi~ (left) and t + t~ (right); the regions below the 
contours are excluded. The black contour ("Super-K up-/i") is the Super-Kamiokande limit 
to 3<r from up-going muons, the orange band is the PAMELA-preferred region, and the red 
ellipses are the Fermi-preferred region; these three are given by Ref. [29J. The dashed green 
line ("AMANDA-II diff.") is the constraint to 90% confidence from the AMANDA-II limit 
on the isotropic diffuse flux of z/ M , and the solid green line ("IceCube case, diff.") is the 
constraint to 90% confidence from the projected IceCube limit on the isotropic diffuse flux 
using cascade events. The dashed blue line ("DeepCore tr. 5yr") is the constraint to 2a 
from IceCube+DeepCore for track-like events after five years of running, and the solid 
blue lines are the constraints to 2a for all-flavor cascade events after one year ("DeepCore 
case, lyr" ) and three years ( "DeepCore case. 3yr" ) of running. 

3 Results 

The results for dark matter decays are shown in Fig. [2] the regions below the contours are 
excluded. The black contour ("Super-K up-//") is the Super-Kamiokande limit to 3a from 
up-going muons discussed in the Introduction. The orange band is the preferred region to fit 
the PAMELA anomaly, and the red ellipses are the preferred region to fit the Fermi e 
anomaly. These three regions are given by Ref. j29] up to mass 30 TeV and lifetime 10 27 s. 

The dashed green line ("AMANDA-II diff.") is the constraint to 90% confidence from the 
AMANDA-II limit on the isotropic diffuse flux of z/ M , and the solid green line ("IceCube case, 
diff.") is the constraint to 90% confidence from the projected IceCube limit on the isotropic 
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diffuse ffux using cascade events. Since the ffux from cosmological dark matter is weak to 
begin with, contributing only 1% of the total flux over the 2-7T facing the galactic center, we 
see that it is quickly overwhelmed by the atmospheric flux below ~ 40 TeV. Nonetheless, 
due to the lower background of atmospheric v e compared to at these high energies, using 
cascade events improves the constraint by a factor of ~ 5. 

The dashed blue line ( "DeepCore tr. 5yr" ) is the constraint to 2a from IceCube+DeepCore 
for track-like events after five years of running, and the solid blue lines are the constraints 
to 2cr for cascade events after one year ("DeepCore case, lyr") and three years ("DeepCore 
case. 3yr") of running. We see that for the final state, while track-like events can only 

reduce the available Fermi-preferred parameter space in five years, cascade events can rule 
it out altogether in only three years. Similarly for the t + t~ final state, track-like events can 
rule out the parameter space in less than five years, but with cascade events it will only take 
one year. Note the weakening of the constraints below m = 250 GeV, where the energy per 
final-state particle is less than 125 GeV. This is caused by the energy cutoff at 40 GeV. 

The results for annihilation are shown in Fig. [3j the regions above the contours are 
excluded. The plots show the same constraints as for decay, except that no isotropic limits 
are shown because they are weaker than the Super-Kamiokande limit by a factor ~ 10 5 due 
to the low density of dark matter on cosmological scales. As with decays, cascade events 
greatly accelerate the development of a useful constraint. For the final state the region 

by the Fermi data can be eliminated in only one year. 

The exclusion plot for annihilations to r + r _ is shown only for completeness, as the 
Fermi-preferred region has already been eliminated by the Super-Kamiokande observation of 
upward-going muons. However, it may provide a useful generic constraint on all dark matter 
models irrespective of the e anomalies. 

4 Conclusion 

We have shown that, by using cascade events, IceCube+DeepCore can more quickly establish 
constraints on dark matter models that would explain the reported e anomalies, and over 
time establish stronger constraints than from track-like events. Specifically, track-like events 
will be able to significantly constrain the parameter space of decays to , and rule out 
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Figure 3: Constraints for annihilation to fi + fi~ (left) and t + t~ (right); the regions above 
the contours are excluded. The black contour ("Super-K up-/i") is the Super-Kamiokande 
limit to 3a from up-going muons, the orange band is the PAMELA-preferred region, and the 
red ellipses are the Fermi-preferred region; these three are given by Ref. [22] • The dashed 
blue line ( "DeepCore tr. 5yr" ) is the constraint to 2a from IceCube+DeepCore for track- 
like events after five years of running, and the solid blue lines are the constraints to 2a for 
all-flavor cascade events after one year ("DeepCore case, lyr") and three years ("DeepCore 
case. 3yr") of running. 

decays to r + r~ and annihilations to fi + fi~ in less than five years of running. In comparison, 
cascade events can rule out decays to yU + /i~ in only three years, and rule out decays to 
t + t~ and annihilation to fi + fi~~ after only one year. Moreover, these constraints are highly 
robust to the choice of dark matter halo profile and independent of dark matter-nucleon 
cross section. 

In closing, we note two interesting possibilities for future work. First, if the pointing 
accuracy for track-like events at IceCube+DeepCore is established to be less than 10°, the 
signal-to-noise for annihilations may be significantly enhanced by observing a smaller region 
around the galactic center, possibly out-performing cascade searches (albeit with greater 
dependence on the choice of dark matter halo profile). This would strengthen the discovery 
potential for dark matter because the galactic center could be identified as a localized source 
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of excess neutrinos. Second, because the v T atmospheric background is so low at energies 
above 40 GeV and at low zenith angles, if IceCube+DeepCore can demonstrate efficient v T 
discrimination [51)], signal to noise could be increased by a factor ~ 100. This would put 
leptophilic dark matter to a severe test. 

5 Acknowledgments 

The authors would like to thank D. Grant, D. J. Koskinen, and I. Taboada for answering 
questions about DeepCore. 

SKM is supported by World Premier International Research Center Initiative (WPI Ini- 
tiative), MEXT, Japan, and would also like to thank G. Lambard for useful discussions. MRB 
is supported by the U.S. Department of Energy, under Grant No. DE-FG03-92-ER40701. KF 
is supported by the U.S. Department of Energy and MCTP via the University of Michigan. 
DS is supported by NSF Grant No. AST-0507117 and GAANN (D.S.). HM is supported 
in part by World Premier International Research Center Initiative (WPI Initiative), MEXT, 
Japan, in part by the U.S. DOE under Contract No. DE-AC03-76SF00098, and in part by 
the NSF under Grant No. PHY-04-57315. 

References 

[1] J. Dunkley, et al, Astrophys. J. Suppl. 180, 306 (2009) [arXiv:0803.0586vl [astro-ph]]. 

[2] G. Jungman, M. Kamionkowski and K. Griest, Phys. Rept. 267, 195 (1996) [arXiv:hep- 
ph/9506380]. 

[3] G. Bertone, D. Hooper and J. Silk, Phys. Rept. 405, 279 (2005) [arXiv:hep-ph/0404175]. 

[4] M. Srednicki, K. A. Olive and J. Silk, Nucl. Phys. B 279, 804 (1987). 

[5] K. Freese, Phys. Lett. B 167, 295 (1986); L. M. Krauss, M. Srednicki and F. Wilczek, 
Phys. Rev. D 33, 2079 (1986). 

[6] J. R. Ellis, R. A. Flores, K. Freese, S. Ritz, D. Seckel and J. Silk, Phys. Lett. B 214, 403 
(1988); M. S. Turner and F. Wilczek, Phys. Rev. D 42, 1001 (1990); M. Kamionkowski 
and M. S. Turner, Phys. Rev. D 43, 1774 (1991); J. Silk and M. Srednicki, Phys. Rev. 



11 



Lett. 53, 624 (1984); L. Bergstrom, P. Ullio and J. H. Buckley, Astropart. Phys. 9, 137 
(1998) [arXiv:astro-ph/9712318]. 

[7] N. W. Evans, F. Ferrer and S. Sarkar, Phys. Rev. D 69, 123501 (2004) [arXiv:astro- 
ph/0311145]; L. Bergstrom and D. Hooper, Phys. Rev. D 73, 063510 (2006) [arXiv:hep- 
ph/0512317]. 

[8] P. Gondolo, Phys. Lett. B 295, 104 (1992). 

[9] O. Adriani et al. [PAMELA Collaboration], Nature 458, 607 (2009) [arXiv:0810.4995 
[astro- ph]]. 

[10] J. Chang et al. [ATIC Collaboration], Nature 456, 362 (2008). 
[11] S. Torii et al, arXiv:0809.0760 [astro-ph]. 

[12] A. A. Abdo et al. [The Fermi LAT Collaboration], Phys. Rev. Lett. 102, 181101 (2009) 
[arXiv:0905.0025 [astro-ph.HE]]. 

[13] D. Hooper, P. Blasi, P. D. Serpico, JCAP 0901, 025 (2009) [arXiv:0810.1527 [astro-ph]]; 
S. Profumo, arXiv:0812.4457 [astro-ph]. 

[14] C. R. Chen, F. Takahashi and T. T. Yanagida, Phys. Lett. B 671, 71 (2009) 
[arXiv:0809.0792 [hep-ph]]; Phys. Lett. B 673, 255 (2009) [arXiv:0811.0477 [hep-ph]]. 
C. R. Chen, M. M. Nojiri, F. Takahashi and T. T. Yanagida, Prog. Theor. Phys. 122 
(2009), 553 [arXiv:0811.3357 [astro-ph]]. 

[15] C. R. Chen and F. Takahashi, JCAP 0902, 004 (2009) [arXiv:0810.4110 [hep-ph]]. 

[16] P. F. Yin, Q. Yuan, J. Liu, J. Zhang, X. J. Bi, S. H. Zhu and X. Zhang, Phys. Rev. D 
79, 023512 (2009) [arXiv:0811.0176 [hep-ph]]. 

[17] K. Ishiwata, S. Matsumoto and T. Moroi, Phys. Lett. B 675, 446 (2009) 
[arXiv:081 1.0250 [hep-ph]]. 

[18] S. Shirai, F. Takahashi and T. T. Yanagida, Phys. Lett. B 675 , 73 (2009) 
[arXiv:0902.4770 [hep-ph]]; arXiv:0902.4770 [hep-ph]. 



12 



[19] A. Ibarra, A. Ringwald, D. Tran and C. Weniger, JCAP 0908, 017 (2009) 
[arXiv:0903.3625 [hep-ph]]; see also A. Ibarra, A. Ringwald and C. Weniger, JCAP 
0901, 003 (2009) [arXiv:0809.3196 [hep-ph]]. 

[20] M. Ibe, H. Murayama, S. Shirai and T. T. Yanagida, JHEP 0911, 120 (2009) 
[arXiv:0908.3530 [hep-ph]]. 

[21] V. Barger, W. Y. Keun, D. Marfatia and G. Shaughnessy, Phys. Lett. B 672, 141 (2009) 
[arXiv:0809.0162 [hep-ph]]. 

[22] I. Cholis, L. Goodenough, D. Hooper, M. Simet and N. Weiner, Phys. Rev. D 80, 123511 
(2009) [arXiv:0809.1683 [hep-ph]]. 

[23] L. Bergstrom, J. Edsjo and G. Zaharijas, Phys. Rev. Lett. 103, 031103 (2009) 
[arXiv:0905.0333 [astro-ph.HE]]. 

[24] A. Sommerfeld, Annalen der Physik 403, 257 (1931); M. Cirelli, A. Strumia, 
arXiv:0808.3867 [astro-ph]; M. Cirelli, M. Kadastik, M. Raidal and A. Strumia, Nucl. 
Phys. B 813, 1 (2009) [arXiv:0809.2409 [hep-ph]]. 

[25] R. Harnik and G. D. Kribs, Phys. Rev. D 79, 095007 (2009) [arXiv:0810.5557 [hep-ph]]; 
A. E. Nelson and C. Spitzer, arXiv:0810.5167 [hep-ph]; I. Cholis, D. P. Finkbeiner, 
L. Goodenough and N. Weiner, JCAP 0912, 007 (2009) [arXiv:0810.5344 [astro-ph]]; 
K. M. Zurek, Phys. Rev. D 79, 115002 (2009) [arXiv:0811.4429v3 [hep-ph]]; P. J. Fox 
and E. Poppitz, Phys. Rev. D 79, 083528 (2009) [arXiv:0811.0399 [hep-ph]]; I. Cholis, 
G. Dobler, D. P. Finkbeiner, L. Goodenough and N. Weiner, Phys. Rev. D 80, 123518 
(2009) [arXiv:0811.3641 [astro-ph]]; S. C. Park and J. Shu, Phys. Rev. D 79, 091702(R) 
(2009) [arXiv: 090 1.0720 [hep-ph]]; D. Hooper and K. M. Zurek, Phys. Rev. D 79, 103529 
(2009) [arXiv:0902.0593 [hep-ph]]; D. Hooper and T. M. P. Tait, Phys. Rev. D 80, 
055028 (2009) [arXiv:0906.0362 [hep-ph]]. 

[26] N. Arkani-Hamed, D. P. Finkbeiner, T. R. Slatyer and N. Weiner, Phys. Rev. D 79, 
015014 (2009) [arXiv:0810.0713 [hep-ph]]. 

[27] G. Bertone, M. Cirelli, A. Strumia and M. Taoso, JCAP 0903, 009 (2009) 
[arXiv:0811.3744 [astro-ph]]. 



13 



[28] S. Palomares-Ruiz, Phys. Lett. B 665, 50 (2008) [arXiv:0712.1937 [astro-ph]]; J. Hisano, 
M. Kawasaki, K. Kohri, K. Nakayama, Phys. Rev. D 79, 043516 (2009) [arXiv:0812.0219 
[hep-ph]]. 



[29 
[30 
[31 
[32 
[33 

[34 

[35 



[36 



[37 



[39 



[40 



[41 



P. Meade, M. Papucci, A. Strumia and T. Volansky, arXiv:0905.0480 [hep-ph]. 

A. Ibarra, D. Tran, Phys. Rev. Lett. 100, 061301 (2008) [arXiv:0709.4593 [astro-ph]]. 

C. -R. Chen, S. K. Mandal, F. Takahashi, arXiv:0910.2639 [hep-ph]. 

O. Adriani et al, Phys. Rev. Lett. 102, 051101 (2009) [arXiv:0810.4994 [astro-ph]]. 

D. Feldman, Z. Liu and P. Nath, Phys. Rev. D 79, 063509 (2009) [arXiv:0810.5762 
[hep-ph]]. 

M. Ibe, H. Murayama and T. T. Yanagida, Phys. Rev. D 79, 095009 (2009) 
[arXiv:0812.0072 [hep-ph]]. 

A. Bhadra and R. K. Dey, Mon. Not. R. Astron. Soc. 395, 1371 (2009) [arXiv:0812.1845 
[astro-ph]]; W. Bednarek, Phys. Rev. D 79, 123010 (2009) [arXiv:0906.0084 [astro- 
ph.HE]]. 

R. Gandhi, C. Quigg, M. H. Reno and I. Sarcevic, Astropart. Phys. 5, 81 (1996) 
[arXiv:hep-ph/95 12364]. 

R. Gandhi, C. Quigg, M. H. Reno and I. Sarcevic, Phys. Rev. D 58, 093009 (1998) 
[arXiv:hep-ph/9807264]. 

J. F. Beacom and J. Candia, JCAP 0411, 009 (2004) [arXiv:hep-ph/0409046]. 

J. F. Beacom, N. F. Bell and G. D. Mack, Phys. Rev. Lett. 99, 231301 (2007) 
[arXiv:astro-ph/0608090]. 

H. Yuksel, S. Horiuchi, J. F. Beacom and S. Ando, Phys. Rev. D 76, 123506 (2007) 
[arXiv:0707.0196 [astro-ph]]. 

A. Gould, Astrophys. J. 388, 338 (1992); G. Jungman, M. Kamionkowski and K. Griest, 
Phys. Rep. 267, 195 (1996); D. Hooper and J. Silk, New Journ. of Phys. 6, 23 (2004). 



14 



[42] M. Kowalski, JCAP 0505, 010 (2005) [arXiv:astro-ph/0505506]. 

[43] L. Gerhardt [IceCube Collaboration], Proc. of the 30th ICRC (2007) [arXiv:0711.0353 
[astro-ph]]. 

[44] J. Kiryluk [IceCube Collaboration], Proc. of the 31st ICRC (2009) [arXiv:0909.0989 
[astro-ph.HE]]. 

[45] S. R. Klein [IceCube Collaboration], IEEE Trans. Nulc. Sci. 56, 1141 (2009) 
[arXiv:0807.0034 [astro-ph]]. 

[46] E. Resconi [IceCube Collaboration], Nucl. Instrum. Meth. A 602, 7 (2009) 
[arXiv:0807.3891 [astro-ph]]. 

[47] D. Spolyar, M. Buckley, K. Freese, D. Hooper and H. Murayama, [arXiv:0905.4764 
[astro-ph.CO]]. 

[48] M. R. Buckley, K. Freese, D. Hooper, D. Spolyar and H. Murayama, [arXiv:0907.2385 
[astro-ph.HE]]. 

[49] E. Middell, J. McCartin and M. D'Agostino [IceCube Collaboration], Proc. of the 31st 
ICRC (2009). 

[50] D. Grant, D. J. Koskinen and C. Rott [IceCube Collaboration], Proc. of the 31st ICRC 
(2009). 

[51] M. Honda, T. Kajita, K. Kasahara, S. Midorikawa and T. Sanuki, Phys. Rev. D 75, 
043006 (2007) [arXiv:astro-ph/0611418]; J. G. Learned and K. Mannheim, Ann. Rev. 
Nucl. Part. Sci. 50, 679 (2000). 

[52] F.-F. Lee and G.-L. Lin, Astropart. Phys. 25, 64 (2006). 

[53] R. Catena and P. Ullio, arXiv:0907.0018 [astro-ph.CO]. 

[54] R. Abbasi et al. [IceCube Collaboration], Phys. Rev. D 79, 102005 (2009) 
[arXiv:0902.0675 [astro-ph.HE]]. 



15 



[55] J. F. Navarro, C. S. Frenk and S. D. M. White, Astrophys. J. 490, 493 (1997) 
[arXiv:astro-ph/96 11107]. 

[56] T. Sjostrand, S. Mrenna and P. Skands, JHEP 0605, 026 (2006) [arXiv:hep- 
ph/0603175]. 

[57] A. B. McDonald, C. Spiering, S. Schoenert, E. T. Kearns and T. Kajita, Rev. Sci. 
Instrum. 75, 293 (2004) [arXiv:astro-ph/03 11343]. 



16 



